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Abstract 

Consider a three-dimensional fluid in a rectangular tank, bounded by a flat 
bottom, vertical walls and a free surface evolving under the influence of gravity. 
We prove that one can estimate its energy by looking only at the motion of 
the points of contact between the free surface and the vertical walls. The 
proof relies on the multiplier technique, the Craig-Sulem-Zakharov formulation 
of the water-wave problem, a Pohozaev identity for the Dirichlet to Neumann 
operator, previous results about the Cauchy problem and computations inspired 
by the analysis done by Benjamin and Olver of the conservation laws for water 
waves. 


1 Introduction 


Consider surface waves over an incompressible liquid, evolving under the influence 
of gravity, in the case where the fluid is located inside a fixed rectangular tank of 
the form M = Q x [—h, -boo) where Q = [0, Li] x [0, L 2 ] and h is a positive constant. 
At time t, the fluid domain U(t) is given by 

U(t) = { (x, y) : X gQ, -h<y < r]{t, x) } , (1) 

where x = (xi,X 2 ) (resp. y) is the horizontal (resp. vertical) space variable. The 
equations which dictate the motion are the incompressible Euler equations with 
free surface. This is a system of two nonlinear equations: the incompressible Euler 
equation for the velocity potential iji: 17 —)• M (so that the velocity is u = 
and a kinematic equation for y which states that the free surface moves with the 
fluid. The energy, which is the sum of the potential energy and the kinetic energy, 
is conserved: 

^ = 0 with n = ^[ y‘^{t,x)dx+ ^ [[ \Vx,y4'{t,x,y)\^ dxdy, (2) 

2, Jq 2 

where g is the acceleration of gravity. This paper is devoted to the analysis of the 
following question: is-it possible to estimate the energy of gravity water waves by 
looking only at the motion of some of the curves of contact between the free surface 
and the vertical walls? From the point of view of control theory, this is the question 
of boundary observability of gravity water waves. 

This work is partly supported by the grant “ANAE” ANR-13-BS01-0010-03. 
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1.1 The water-wave equations 


Hereafter we use the notations V = (i9a;i, '^x,y = i^x,dy) and = A^ + dy. 

As already mention, the velocity field u: 0 —>■ is given by u = ^x,y4> for some 

velocity potential (/>: 11 —M satisfying 

Ax,y4> = + - |Va:, 2 y(/>|^ + P + 52/= 0, (3) 

where P is the pressure term. The water-wave equations are then given by several 
boundary conditions. Firstly, the velocity satishes the solid wall boundary condition, 
so that dn4> = 0 on the boundary of the tank where dn denotes the outward 


normal derivative. This implies that 

dxi4> = 0 for xi = 0 or xi = Li, (4) 

5x2 </> = 0 X 2 = 0 or X 2 = L 2 , (5) 

dy4> = 0 for 2 / = —h. (6) 

The problem is then given by two boundary conditions on the free surface {y = i]}- 

dtij = s/TTW^dn(j)\y=ri = dy4>{t,x,r]) - Vr]{t,x) • V4>{t,x,r]), (7) 

^\y=v ~ (^) 


The equation (7) implies that the free surface moves with the fluid. The condition 
P{t, x,r]) = 0 is a dynamic condition that expresses a balance of forces across the 
free surface. 


We also assume that the free surface must intersect the vertical walls orthogonally: 


dxiV = 0 for xi = 0 or xi = Ti, 
dx 2 V = 0 foi' X 2 = 0 or X 2 = L 2 . 


(9) 


It is proved in [2] that (9) always holds for smooth enough solutions (we elaborate 
on that fact in Section §2.4). Eventually we assume that 


r] > 




= 0 . 


( 10 ) 


One can always assume that the mean value vanishes since it is a conserved quantity. 


1.2 Boundary observability of water waves 

There are very few articles about the possible applications of control theory to 
the study of water waves. The first results are due to Reid and Russell [20] and 
Reid [18, 19] who studied the linearized equations at the origin. Alazard, Baldi 
and Han-Kwan initiated in [1] the study of the control of the nonlinear equations. 
The analysis in [1] and the one in this paper rely on completely different tools. 
This is because of the following differences: the article [1] addresses the internal 
control problem for the two-dimensional equations with surface tension, while here 
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we consider the boundary observability for the 2D and 3D problems without surface 
tension. In [1], the analysis used in an essential way the infinite speed of propagation 
of gravity-capillary waves, which is the property that an harmonic with frequency |^| 
travels at a speed proportional to By contrast, for gravity waves, an harmonic 

with frequency |^| travels at a speed proportional to Loosely speaking, in 

this paper we study the observability problem for low and medium frequency gravity 
waves, while in [1] the controllability of gravity-capillary waves is deduced from an 
observation of high frequencies. With regards to the proofs, the paper [1] relies on 
microlocal analysis while here the proof will rely on the study of global quantities. 
More precisely, in [1] the analysis combines Ingham type inequalities, paradifferential 
calculus and other techniques used to study quasi-linear problems. By contrast, in 
this paper, the analysis will be based on the multiplier method and on various exact 
identities which are studied as conservations laws. In this direction, let us mention 
the recent article by Biccari [8] which introduced the use of the multiplier method to 
analyze the interior controllability problem for the fractional Schrodinger equation 
idtu+{—A)^u = 0 with s > 1/2 in a bounded domain Q with Dirichlet boundary 
condition^. 


We do not assume that the reader is familiar with control theory and begin by 
recalling what is the multiplier method in the simplest case. Consider the one 
dimensional linear wave equation with Dirichlet boundary condition: 

(9/tt — = 0, u{t,0) = u{t,l) = 0. (11) 


Multiply the equation by xdxU and integrate by parts, to obtain 

^ [ {dxu{t,l)f dt = [ {dtu){xdxu) dx + 7 ; [[ [(dtu)'^ + (d^uf] dxdt (12) 
^ Jo Jo 0 ^JJs 


where S = (0,r) x (0,1). Since 

fi 

{dtu){xdxu) dx 


f 


<S:=5 


[{dtuf + {dxuf] dx, 
by using the conservation of energy {dS/dt = 0), we deduce 

{dxu{t,l))‘^ dt> {T-2) f [{dtuf + {dxuf]{0,x)dx. 
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Jo 


(13) 


(14) 


This inequality implies that, for T > 2, one can bound the energy by means of an 
observation at the boundary. 


There are inequalities analogous to (14) which hold in multi-dimensional domains. In 
this direction, let us also mention the celebrated result by Bardos-Lebeau-Rauch [6] 
which is devoted to the boundary observability of solutions of second-order hyper¬ 
bolic equations in the general case when the region of observability meets every ray 
of geometric optics. 

^ One cannot apply this result to study the observability of gravity water waves. Indeed, i) we 
are concerned with a nonlinear problem and ii) even the linearized problem involves a different 
setting since it can be written under the form idtu -\- { — A)‘‘u = 0 with s = 1/4 (and hence the 
assumption s > 1/2 does not hold). The key point is that, if s < 1/2, then high frequency waves 
propagate at a speed which goes to 0 when |^| goes to -too. 
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The inequality (14) is a central result in the study of the control of the wave 
equation—extensions and applications of (14) are discussed in the first part of the 
SIAM Review article by Lions [16]; see also [11, 23, 17, 22]. In particular, (an ex¬ 
tension of) (14) allows to prove that in multi-dimensional domains, one can drive a 
solution of the wave equation from the rest position to a desired state, in finite time, 
by acting only on part of the boundary. Our intention in this paper is to start the 
study the analogous problem for the water-wave equations. This is a very natural 
question since it corresponds to the wavemaker problem: the aim is to determine 
which waves can be produced by, say, the motion of a plate immersed in a fluid. 


In this paper, our goal is to obtain a boundary observability result similar to (14) for 
gravity waves. By contrast with the wave equation, the water-wave equations are 
fully nonlinear and contain nonlocal terms. However, we will prove that a similar 
observability inequality holds. 


Definition 1.1. Introduce = 4>{t, x,r](t, x)) (evaluation of the potential at 

the free surface) and 

0 ;= -r]dti>- l??^- 
Set 

^(T) = L 1 L 2 j J Q{t,Li,X2)dx2 + ^ j Q{t,xi,L2) dxi^ dt. 

(It is convenient to use the factor L 1 L 2 to derive some identities.) 


Notice that ^{T) depends only on the boundary values of the unknowns. Our main 
result asserts that 3S{T) is larger than the energy % (cf (2)) when T is large enough. 

Theorem 1.2. Let j3 > 1/2 and x be either a C'“(Q) function or y = 1. There 
exist three positive constants Kq, k, c such that, for any N in N, the following result 
holds. Assume that the initial data rjo and ijjQ are of the form 


Xi^) ^ ^ CLnm 


COS 


nx\ \ 

vr— \ cos 


vr 


mx2 \ 


with jarimi < cN 


(15) 


and such that rjo satisfies (10). Then there exists a unique smooth solution (ry, (/>) of 
the water-wave equations (3)-(10) such that ry(0,x) = ryo(x) and '4){0,x) = fio{x). 
This solution exists on the time interval [0,T7v] with Tjq = Kq KqN^^ and satisfies 


^{Tn) > n. 


If, in addition, rj and 4> are independent of X 2 (this means that we consider a two- 
dimensional wave), then one has the stronger conclusion 

cTm 

Li / Q{t,Li,0)dt >71. 

Jo 


We refer to Section 2 for an explanation of what we call a smooth solution of the 
water-wave equations. Several other remarks are in order. 
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(i) One can consider more general initial data, see Remark 5.3. 


(ii) The second point to be made is a clarification of how one passes from a true 
observation at the boundary to a control of ^{T) or Q{t, Li,0) dt. Let us 
explain why these two quantities can be expressed as quantities depending only the 
restrictions to dQ of r] and dtr] (notice that ry, dtrj are quantities which can measured 
by a camera by contrast with ^). To see this, consider firstly the case of a two- 
dimensional wave and set m{t) = Then one has 

0(f, Ml, 0) = i , 

as can be verified by a simple calculation (see Remark 4.6). For a three-dimensional 
wave, this is not so simple. However, one can determine ^|J from (ry, dtr]) by using the 
equation (7) and considering the Neumann to Dirichlet operator (as in Appendix 
A.3 in Lannes’ book [15]). 


(in) Theorem 1.2 will be deduced from an observability estimate which holds without 
smallness assumptions. Namely, we will prove that for all smooth enough solution 
defined on the time interval [0,T], one has (cf (51)) 


^(T) > 



5 +2d 


BTj n 



max{Li,L2} 




A 


where 

B := sup ||V7y(f)||^oc, A := sup ||VV'(f)||^ 2 . 
te[o,T] te[o,r] 

The assumption (15) will be used only to prove that one can bound A in terms of 
y/n on large time intervals. Then, when T is large enough, the right-hand side of 
the above inequality is larger than %. 


(iv) Theorem 1.2 gives an observability inequality in time for /3 > 1/2. The 
important think to note is that, except for the fact that we preclude j3 = 1/2, this 
result is sharp. Indeed, an harmonic with frequency N travels at a speed proportional 
to 1 /y/N (the dispersion relationship for the linearized equations is oj(k)‘^ = g\k\). 
So it might take a time y/N to reach the boundary and hence for the observation to 
be possible. Now we can explain why we need a smallness assumption on the initial 
data, namely the assumption |a„m| < cN~'^. This is because one needs such a 
smallness assumption to guarantee that the solutions exist on large time intervals of 
size (the link between n and /3 will be made through an interpolation argument 
in the proof). 


(v) Another point should be added: as indicated on Figure 1 below, we are making 
an observation on part of the boundary only (satisfying the so-called geometric 
control condition). Indeed, to determine Q{t, Li,0) dt it is sufficient to look 
at the motion of the point of contact ^(t) = (Li, 0,7y(f, Li, 0)) between the free 
surface and only one wall, namely {xi = Li}. Similarly, for a three-dimensional 
wave, to determine AS{T), it is sufficient to observe the motion of the curves of 
contact between the free surface and the walls {xi = Li} and {x 2 = L 2 }: 

^i{t) = {(Li,X2,y), X2 G [0,7.2], y = ry(t,Ti,X2)}, 

'^ 2 {t) = {(xi,L 2 ,y), xi G [0,Li], y = 7y(f,xi, 7 . 2 )}- 
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Figure 1: Three-dimensional and two-dimensional waves in a rectangular tank 

1.3 Strategy of the proof and main identity 

In the rest of this introduction, for the sake of simplicity, we consider two-dimensional 
waves (these are solutions independent of X 2 , and we write simply r]{t,xi) and 
xi,y)). We refer to the next sections for the corresponding statements for three- 
dimensional waves. 

Theorem 1.2 will be deduced from an exact identity, similar to (12), where, quite 
surprisingly, the terms coming from the nonlinear part of the equations have a very 
simple and compact form. 

Theorem 1.3. Consider a smooth enough solution of the water waves equations 
(3)-(10) defined on the time interval [0, T] and independent of X 2 - Introduee 


m{t) = r]{t,Li). 

Then Q{t,Li) = ^[gm{t)‘^ — m(t)m'(f)^] and 






where Ti is given by (2) and J f dx\ 


t=T 

t=o 


stands for f f{T, x) dx — f /(O, x) dx. 
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Theorem 1.3 is proved in Section 4. The proof uses the Zakharov’s formulation 
of the water-wave problem as a Hamiltonian system (see [24]) and the observation 
by Craig and Sulem [12] that the equations and the hamiltonian are most naturally 
expressed in terms of the Dirichlet to Neumann operator G{r]). The main ingredients 
of the proof of Theorem 1.3 are then: i) a Pohozaev identity for the Dirichlet 
to Neumann operator (that is a computation of f {G{ri)'il))xdx'>p dx) which shows 
that the contributions due to the boundary conditions are positive and ii) some 
computations inspired by the analysis of Benjamin and Olver [7] of the conservation 
laws for water waves. Let us mention that, in the appendix, we give another proof 
of (16) which exploits the hamiltonian structure of the water-wave equations. The 
idea is to compute 


A := 


[0,T]x[0,Li] 


[(dt?])(xdx'tp) - (dt'ip)(xdxrj)] dxdt. 


We compute A in two different ways, and the wanted identity (16) is obtained by 
comparing the two results. The first computation is simply an integration by parts. 
The second computation relies on the fact that the equation is hamiltonian (see 
Zakharov [24], Craig-Sulem [12]), which means that 


Therefore 





m 

6r] ' 


A = 



with 



, m 5% 

+ XTjx — 

dip dr] 


dx, 


Then we compute a by writing 


a = lim - + exipx) — T~i{ri, ■;/’)] -|- lim - \H{ri -|- £xrix,ip) — T-L{rj, ?/i)] , 

£^0 £ £->0 £ 


and using a shape derivative formula due to Lannes [14]. This allows to avoid the use 
of the Pohozaev identity. However, the proof still requires to prove some identities 
for quantities which are analyzed as conservation laws. 


We now indicate how to infer from Theorem 1.3 a boundary observability result. 
This will require to make an additional assumption. Indeed, in sharp contrast with 
the example of the wave equation (11) discussed above, even for the linearized water- 
wave equations at the origin, there is no observability inequality in hnite time which 
applies to all initial data (see Reid and Russell [20] or Biccari [8] for related results). 
However, we will prove a boundary observability result, for the full nonlinear prob¬ 
lem, for initial data whose Fourier transforms are compactly supported. To do so, 
we begin with the following corollary of Theorem 1.3. 


Corollary 1.4. Consider a smooth enough solution of the water waves equations 
(3)-(10) defined on the time interval [0,T]. Assume that there exists two positive 
constants A, B such that 


B < 


2 

7’ 


T > 


4 

2 - 7R 



5Li 

7^' 


A 


(17) 
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and 


(18) 


sup \\da;V{t)\\Lo.mLi])<B, sup \\d^i;{t)\\^ 2 noL^]^<AVm. 

t€[0,T] tg[0,T] 


Then 


Li / eit,Li)dt>n 
Jo 


T pT rm(t) 

+ ^ / / {dy(t)f{t,Li,y)dydt 

^ Jo J-h 

h 

+ / / idx(pf{t,x,-h)dxdt. 

16 Jo Jo 


(19) 


Proof. The proof is straightforward: directly from (18), the definition of TJ, the 
Cauchy-Schwarz inequality and the assumption rjdx = 0, we get 


rL\ j-Lx / j-Lx \ 

/ yf) dx = / V {'4’ ~ f) dx\ dx 


< Li M\lI 


1^2 < {2L,/^)An, 


4 Jo J Jn{t) 


idx'n){dxj))idy(j)) dxdydt 


7B 

<-Tn, 


fL, 

/ xyd^ 
Jo 


dx 


< Li M^2 \\dx4^\\L2 < {2L,/^)An, 


where we used the Poincare inequality in the first line. Now, since rj > —/i/2 by 
assumption, using the identity (16) we see that the left hand side of (19) is larger 
than the sum of the last two terms in the right hand side of (19) and of 


1 5Li , 7B^ , 

-T-—A--T)n. 


which, by assumption (17), is larger than TJ. 


□ 


Remark 1.5. i) One controls not only the energy TJ but also the traces of the 
velocities on the wall {x = Li} and on the bottom. 


ii) {Unique continuation) If m = 0 then (19) implies that TJ = 0 and hence y = 0 
and 0 = 0. 


Hi) The assumption \dx'r]{t,x)\ < 2/7 is physically realistic. Indeed, one expects a 
steepness-induced blow-up in finite time when the wave height is large compared to 
the wavelength (there are no mathematical proof of this claim but it is known that 
blow-up occurs in finite time for some large enough initial data, see [10]). Moreover, 
the threshold 2/7 is in good agreement with the sharp experimental studies of steep 
irregular wave events reported in [13]. 


To conclude this introduction, we explain how we deduce Theorem 1.2 from Corol¬ 
lary 1.4. The reason why we need the assumption for 0 in (18) is the following: we 
do not have an estimate analogous to (13) for the term xydx'f’dx which appears 
in the right-hand side of (16). Indeed, the fact that y/TH is larger than the L‘^{Q.{t))- 
norm of '7/x,y4>{t) gives only a bound of the Hx^^-norm of the trace if. Hence y/TL 










does not control the L^-norm of dx'4’ and this why we need the assumption for ijj in 
(18). Now, notice that this assumption for dxip holds at t = 0 with A = Ky/N if the 
Fourier transform of V'(O) is supported in [—A^] (as in (15)). Since the equations 
are nonlinear, we cannot assume that the Fourier transform of the solution will be 
supported in [—A^, A^] for all time. However, we shall see that, for small data, one 
can propagate the estimate ||<9a;V'(0)||2^2 < KyfN on large time intervals 
and hence deduce and observability result. To do so, we shall combine an interpo¬ 
lation argument and the fact that the Cauchy problem is well-posed on (with s 
large enough) on large time intervals for small initial data. 


Organization of the paper. 

In Section 2 we recall various results about the well-posedness of the Cauchy problem. 
This will allow us to clarify what we call a smooth enough solution. In Section 3, 
we prove a Pohozaev identity for the Dirichlet to Neumann operator and we use 
this identity in Section 4 to prove Theorem 1.3. Then we prove Theorem 1.2 in 
Section 5. In the appendix we give a proof of Theorem 1.3 which exploits the 
hamiltonian structure of the equations. 
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2 About the Cauchy problem in a rectangular tank 


We recall here various results about the well-posedness of the Cauchy problem. 


2.1 The Craig—Snlem—Zakharov system 

To study the Cauchy problem for the water-wave equations, there are at least two 
difficulties. Firstly, one has to reduce the analysis to a time independent domain. 
In the case where the space variable x belongs to the torus instead of Q, several 
approaches have been used. In this paper, we use the Eulerian approach, following 
Zakharov [24], Craig-Sulem [12], Lannes [14, 15] and Alazard-Burq-Zuily [4, 2]. In 
this approach, one works with the trace of the potential cj) at the free surface and 
the Dirichlet to Neumann operator G{r]). 

Hereafter we consider either the 2D problem or the 3D problem (see figure 1) and 
denote by d G {1, 2} the dimension of the free surface. For a three-dimensional fluid 
one has d = 2 and we use the notation V = (d^^,d^jj) (and we also use the notation 
V = dx when d = 1). 
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We set 


= (j){t,x,r]{t,x)), 

and introduce the Dirichlet to Neumann operator G{r]) defined by 

G(??)V’ = y/l + |Vr/|2 dn4)\y=r^ = {dy4> - Vr/ • V^)\y=n 

(see the next paragraph for a precise definition of G{r])). Then the Craig-Sulem- 
Zakharov formulation of the water-wave equations reads 

dtrj = G{r])^, 

+ + — I „ (G(r?)V> + Vr?-VV’)^ = 0- 

^ 2(l + |Vr/| ) 

Moreover, the energy % satisfies 

^ = 1 J ri^dx + ]^ JJ dydx = ^ j [gr]"^ + V’G'(r?)V'] dx. 

We refer the reader to the original article by Craig-Sulem [12] for the derivation 
of this system from the water-wave equations (3)-(8) as well as to [3] for the proof 
that, conversely, given a solution of (20), one can define a solution of (3)-(8). 


2.2 Definition of the Dirichlet to Nenmann operator 

We have already written the water waves equations under the form of the Craig- 
Sulem-Zakharov system (20). We give here the precise definition of the Dirichlet to 
Neumann operator that is used in that system in the context of periodic functions. 

For s E [0,-|-oo), we denote by the Sobolev space of periodic functions: 

= |w = X] ^rien{x ); ^ (1 + \n\'^y \anf < +oo|. 

where en{x) = exp if d = 1, e„(x) = exp “I—Z~)) d = 2 

(we fix the periods equal to 2Li, 2 L 2 for reasons that will be clear below). 

Fix h> 0 and consider g E with sq > d/2 + 1 such that t]{x) > —h/2. Set 

D = {{x,y) E X R : —h <y< r/(x)}. (21) 

Since y is Lipschitz, by the usual variational method, one obtains that, for any 
'll; E the problem 

^x,y4^ — 0 in D, (l)\y=ri — V’) dy'P\y=—h — 0, (22) 

has a unique variational solution. Moreover, even if Vx,y4> belongs only to L^(D), 
one can prove that 

G{v)i> = (dyCf) - VT? • V(l))\y=y, 

is well-defined and belongs to More generally, one has the following 

result (see [4, 15]). 
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Proposition 2.1. Let sq > d/2 + 1 and assume that rj is in Then G{r]) 

is a bounded operator from into for any a E [1/2, sq]- 

Later we shall perform various integrations by parts and apply repeatedly the 
Green’s identity. The fact that all the computations are meaningful relies on the 
following regularity result. 

Proposition 2.2. Letd > 1, s > d/2+2 and assume that E 
Then _ 

Proof. We explain how to deduce this result from the results in Lannes’s book [15]. 

It is convenient to flatten the free surface. To do so, consider a diffeomorphism from 
X [—h, 0] to Cl, of the form (a:, z) !->■ (x, p{x, z)). The simplest choice would be to 
set 

p{x,z) = [l + ^r]{x) +z. (23) 

However, since we only assume that s > d/2 + 2 (the analysis is simpler for s > 
d/2 + 5/2), following Lannes [15, §2.2.2], we need to consider a regularized version 
of (23). Introduce 

C(x, 2 ;) = c Y] e“^^^l"l^r)nen(a;), dn = / en{y)ri{y) dy, 

where c and 5 are positive constants (with 5 small enough), chosen so that 

C|^=o = d, \dzC{x,z)\ f x [-h,0]). (24) 

Then set 

p{x,z) = (l +C{x,z) + z. (25) 

Notice that p{x, 0) = p{x) and p{x, —h) = —h. Recall that p > —h/2 by assumption. 
In view of the bound for dzC, in (24), we deduce that f > —3/i/4 and dzP > 1/5, 
which proves that {x,z) i-)- {x,p{x,z)) is a C^-diffeomorphism. The problem thus 
reduces to establishing that the function ip{x,z) = 4 >{x, p{x, z)) satisfies Vx,zT £ 
C'^(T'^ X [—h, 0]). To obtain this result, we view as a function of z with values in 
functional spaces. Since s > d/2 + 2 and since we chose a special change of variables 
where p is given by (25), we are in position to apply Corollary 2.40 in Lannes’ book 
[15] (with s replaced by s — 1/2). We deduce that 

E L2([_/j,o];dd^-^(T''))nddl([-h,0];F^-i(T'')). 


As a result 

VxT^Ll{[-h,t)\,H^-h{T^)) and 9,E L2([-d, 0]; (¥'’*)), 

dz+£Ll{[-h,t)\-W-" 2 {T’^)) and 5,9,+E L2([_/i, Q]; (T'")), 

and hence 

Vx,zT^Cl{[-h,^]-,H^-\T'^)). ( 26 ) 
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Then, using the second order equation satisfied by (/? one can express d^ip in terms 
of {dzip, dzVxP^, and in terms of p and hence, using the standard nonlinear 

estimates in Sobolev spaces, one deduces that 

5,VgC°([-/i,0];F^-2(T'')). (27) 

The wanted result x [—/i,0]) then follows from (26)-(27) and the 

Sobolev embeddings C C^{T<^) and C C'°(T'^). □ 


2.3 The Cauchy problem for periodic functions 

We recall here a well-posedness result for the Cauchy problem for periodic functions. 
There are now quite a lot of papers on this subject and we quote below only two 
results which are related to our problem. The main difficulty in the analysis of 
System (20) is that writing energy estimates on the function (t/, ijj) makes appear an 
apparent loss of half a derivative. A way to circumvent that difficulty is to bound 
the energy not of (?/,V’) but {r],B,V) where 

i? = V = V(l^\y=, = V^I;-BVp. (28) 

If initially po and V’o belong to 77^o+5(T‘^) for some sq > {d + l)/2, then G{po)’ijjo 
belongs to i7®°“2(T'^) (cf Proposition 2.1). On the other hand, is an 

algebra for any sq > {d + l)/2. It thus follows from usual nonlinear estimates in 
Sobolev spaces that 

^ + Vr?o • VV^o ^ ^ ^ 

1 + IVryor 

(29) 

The following result shows that one can propagate the fact that Bq and Vq are in 
for So > d/2 + 1 (and this is the key point to circumvent the apparent loss 
of half of derivative in the study of the Cauchy problem). 

Theorem 2.3 (from Alazard-Burq-Zuily [4]). Let d > 1, sq > (i/2 + 1 and consider 
an initial data (f?o,V’o) such that 

{Vo,^o,Vo,Bo) G X®'' ;= x i7"«+i(T'^) x x (30) 

Then there exists a time T > 0 such that the Cauchy problem for (20) with initial 
data (po^tpo) has a unique solution {^, 1 /) such that {p,'i/),V,B) G C'°([0,T];X^°). 


2.4 Extension to periodic fnnctions 

We recall here from Alazard-Burq-Zuily [2] how to solve the Cauchy problem for 
the water-wave equations in a rectangular tank. Let (i G {1, 2}, Q = [0, Li] if d = 1 
and Q = [0, Li] x [0, L 2 ] if d = 2. One denotes by v the outward unit normal to Q 
{v = (1,0) if xi = Li, = (0, -1) if X2 = 0,...). 
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As recalled in the introduction, the key observation is the following: for smooth 
enough solutions, the angle between the free surface and the vertical boundary of 
the tank is a right angle (see Section 6 in [2]). This means that di^rj = 0 on dQ. 
Now observe that Vtp = (V<i))\y=r] + {dy<i))\y=ri^rj. Since dn4> = 0 on the boundary 
of the tank ^ = Q x [—/i, +oo), we conclude that = 0 on dQ. 

Definition 2.4. Given l<d<2 and a > 3/2, one denotes by Hg{Q) the space 
H^iQ) = {v e : d^v = 0 on dQ}. 


Consider the Cauchy problem for initial data V’o: Q M in H^{Q) for some 
a large enough. Following Boussinesq (see [9, page 37]), the idea is that, in a 
general setting, one can extend these initial data to periodic functions defined for 
X G solve the Cauchy problem for these extended initial data and then obtain 
a solution to the water-wave equations in a canal by considering the restrictions of 
these solutions. 

Definition 2.5. Let u: Q —)• M. If d = 2 , we define u: —)• M as the unique 


extension of v satisfying 

v{x) = v{x) \/x G Q, (31) 

v{ — Xl,X2)=v{xi,X2)=v{xi,—X2) VxGM^, (32) 

v{xi + 2 Li,X 2 ) = v{x) = v{xi,X2 + 2L2) Vx G M^. (33) 

Similarly, when d = 1, u : M —>■ M is defined by 

v{x) = v{x) Vx G Q, (34) 

v{—x) = v{x) Vx G M, (35) 

v{x + 2 Li) = v{x) Vx G M. (36) 


Definition 2.6. Given cr G M, denote by the Sobolev space of those periodic 

functions which are even (satisfying (32)-(33) when d = 2 and (35)-(36) for d= 1). 

Now consider the case d = 1 (to hx notations) and u G Hf (T) with a > d/2+1 = 3/2. 
Then, dxu{x) is and odd which implies that cia;^i(0) = 0. Moreover, one has 
u{Li + e) = u{—Li + e) = u{Li — e) and hence one has also dxu{Li) = 0 (then 
dxu{nLi) = 0 for any n G Z). We have a similar result when d = 2. This proves 
that 

> ^ + 1, Vx G dd,"(T''), v\QeH:{Q). (37) 

Conversely, the following result shows that any function v in H^{Q), with a G 
(3/2, 7/2), is the restriction to Q of a function belonging to ddg (T'^). 

Proposition 2.7 (from Prop. 6.5 in [2]). Let 1 < d <2 and | < a < |. Then the 
map V t-^v is continuous from Hf{Q) = {u G H^{Q) : dyV = 0 on dQ} to 

We are now in position to define G{rj)ijj when rj and belong to some space Hf{Q). 
To do so, let ?7, V’ be as given by Definition 2.5. If 3/2 < a < 7/2, it follows from 
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Proposition 2.7 that ??, ■0 belong to (T'^). If one further assumes that a > d/2 + 1, 
as recalled in §2.2, there exists a unique periodic variational solution to 

=0 in {(x, y) xR : -h < y < 7/(x)}, 

_ _ _ {6o) 

(/{x, y)\y=+x) = i^{x), dy(l)\y=_h = 0. 

Definition 2.8. Consider y and i/ in H^{Q) with d/2 + 1 < a < 7/2. We define 
G{r])'tp (resp. 4>) by taking the restriction to Q (resp. Q): 

4> = fi\n, G{y)'ilj = G{rj)^\Q. (39) 

where It = {{x,y) & Q x R : —h <y< yix)}. 

Proposition 2.9. Assume that G H^{Q) x H^{Q) with a G (d/2 + 2,7/2) 

where d = 1,2. Then 

Vx,yfi £ C'^(^)) dn4> = 0 on dM n dCl, (40) 

and 

G{y)i,eHr\Q). (41) 

Proof. Let us prove (40). The fact that Vx,y(p G (7^(12) follows from Proposition 2.2. 
It remains only to prove that dnfi = 0 on dM n (90. By definition of fi, we have to 
prove that dn4>{x,y) = 0 for any {x,y) G 0 with x G dQ. To simplify notations, 
assume that d = 1 and Li = 1. Then, as already mentioned after the statement 
of Definition 2.6, notice that rj{x) = rj{—x) and y{Li — x) = rj{Li + x) since rj 
is 2Li-periodic and even. Since a > 3/2, one has rf G H^(T) C C^(T)^and one 
deduces that dx^ifif) = 0 = dxrj{l). Similarly fi{x) = 'ijj{—x), fi{Li — x) = fi{Li + x) 
and, by uniqueness of the harmonic extension, satisfies the same symmetries in x. 
Consequently, dx4> is odd in x and, since dxfi is continuous, we infer that 

dxfi{0,y)=0 Vy G [-/i,r/(0)], dxfi{Li,y) = 0 Vy G [-h,y{Li)]. 

This completes the proof of (40). 

It remains to prove (41). It follows from Proposition 2.1 and Proposition 2.7 that 
G{rf) 2 f G 77°^“^(T'^). Therefore, in view of (37), it remains only to prove that 
G{rj)'ip is even (see Definition 2.6). This in turn follows from the definition G{?j)ip = 
dyfi — S/rj • V(/ I y_~ and the symmetries of (f> that we already used. □ 


2.5 The Cauchy problem in a rectangular tank 

We are now in position to define what is a solution of the equations in a rectangular 
tank and to prove the existence and uniqueness of such solutions. 

Definition 2.10. Let d G {1,2} and s > d/2 + 2 and consider 

(+V’)GC°([0,r];77/(Q)x77/(g)). 


14 


Then (r/, tp) is a solution of the Craig-Sulem-Zakharov system (20) with initial data 
(%)V'o) if and only if the extended functions {rj,'tp) satisfies 

' dtrj = G{rf)\f, 

< dt^ + gv + +|v^|2) (g(^)^ + = 0, (42) 

, {rj,fj)\t=o = 


Proposition 2.11. Let d G {1,2} and s € {d/2 + 2,7/2). Consider initial data 
r|o,^po in Hf{Q) and denote by rjo,'ipo the extensions as given by Definition 2.5. If 

+ Po = V^o-5oV%G^^“^(T‘'), (43) 

1 + IVryor 

then there exists T > 0 and a unique solution {r/,ip) G C^{[0,T]; Hf{T'^) x Hp{T‘^)) 
to the Cauchy problem (42) such that 

{B,V) G C°([0,r];77"-5(T^) x 


where 


~ G(i7)V' +Vi/-VV' 

1 + |Vi/|2 ’ 


V = VV’ — BVrj. 


Proof. In view of Theorem 2.3 it remains only to prove that rj and ip are even in 
for 1 < A: < d (satisfying (32)-(33) when d = 2 and (35)-(36) for d = Lj. To do so, 
assume that d = 2 and k = 1. Set r]^{t,x) = r/(t, —xi,X 2 ), ip'^{t,x) = ip{t,—xi,X 2 ). 
Since the equations are invariant by the symmetry xi 1 —>■ —xi, we get that {ri'^,ip'^) 
satisfies the same equations. Moreover, one has r]'^\t=o = Vo, 'P>Ht=o = 'Po since r/o 
and ipo are even in Xfc. By uniqueness, we deduce that = rj and ip'^ = ip, which is 
the desired property. □ 

Remark 2.12. Let us give three cases where the assumptions are satisfied, i) r/o, ipo 
are finite linear combinations of terms of the form (15). 

ii) Consider the case where i/jq = 0. Then ipQ vanishes and hence Bq = 0 and Vq = 0 
so Assumption (43) holds whenever r/o G R®(T'^) (and hence for any r/o G Hf{Q)). 

in) Assume that d =^1 and r]Q,ipo G Hf{Q) for some a G (3,7/2). Proposition 2.7 
implies that r/o and r/o belong to Rg (T'^)^. Since rr > 3, it follows from Proposi¬ 
tion 2.1 below and the usual nonlinear estimates in Sobolev spaces that Bq and Vq are 
in 77°"“^(T'^). So we may apply the assumptions of Corollary 2.11 with s = a — 1/2. 


3 Pohozaev identity 


In this section the time is seen as a parameter and we skip it. As above, we denote 
by d G {1,2} the dimension of the free surface. For a 3D (resp. 2D) fluid one has 
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d = 2 (resp. d = 1) and we use the notation V = (9a:n (resp. V = dx)- Our goal 
is to prove a Pohozaev type identity for G{r])'tp, that is we want to compute 

f {G{ri)'ip){xdx. 

Jq 

Recall from the previous section that G{iri)^ is defined by taking the restrictions to 
Q of a periodic function G{?j)^jJ (see Definition 2.8). Consequently, even if we are 
working in a bounded domain, we are essentially handling periodic functions defined 
on R'^. Recall also that one has dn4> = 0 on n (9D (see Proposition 2.9). 

Proposition 3.1 (Pohozaev identity). Assume that (r?,?/') belongs to H^{Q) x 
H^{Q) for some a > d/2 + 2. Denote by R the solid part of dD: 

R := dM n dD, 

= Q X [—h,+oo)) and denote by n the unit outward normal to dD. Then, 
f {G{r])ip){x-Vip) dx 

Jq 

i^y) dxdy (44) 

+ ^ [ [r] — X ■ Vr/) — 2BG{r])fj] dx, 

2 Jq 

where (j) is given by (39), B = {dy4>)\y=^(^x) ^ x4>)\y=y{x) (Proposition 2.9 

implies that Vx,y4> G (7^(14) and hence all the terms are well-defined). 

Remark 3.2. f) If d = 1 then the second term in the right-hand side of (44) vanishes 
and, since n • = 0 on R, the first one simplifies to 

1 r /'T’X T h 

2 • ndS = ^ {dy(p{Li,y))‘^ dy + - {dx(j){x,-h)f dx. 

Consider now the case d = 2. Then one has also (x, y) • n > 0. Indeed, 

on {xi = Li} one has n = (1,0,0) and (x, y) ■ n = Li, 

on {x 2 = T 2 } one has n = (0,1,0) and (x, y) ■ n = L 2 , 

on {y = —h} one has n = (0,0, —1) and (x, y) ■ n = h, 

and moreover, (x, y) • n = 0 on the two other faces {xi = 0} and {x 2 = 0}. 


ii) Another Pohozaev identity for the fractional Laplacian has been proved by Ros- 
Oton and Serra [21]. Moreover, Biccari [8] deduced from the Pohozaev identity in 
[21] a Pohozaev identity for solutions to fractional Schrodinger equations. Namely, 
it is proved in [21] that, if Q is any domain of R'^, s G (0,1) and u G 
vanishes in R'^ \ Q, then 



—— / u{—Ayudx 
2 Jq 

r(i + s)^ r / u 

2 JaQ Vdist(x,aQ)^ 


2 

(x • n) dS. 


(45) 
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To compare both results the important think to note is that, in the case without 
boundary, one has G(0) = (—However, (44) and (45) involve functions sat¬ 
isfying different boundary conditions. Another essential difference for our purpose 
is that (44) applies in the variable coefficients case where rj ^ 0. Also the proofs of 
(44) and (45) are different. The proof of (44) given below is in fact guided by the 
study of the commutator [xdx, G{r])] in Alazard-Delort [5, Chapter 4]. However, one 
cannot apply the results of [5] because of the boundary conditions on R (and also 
because we consider the case d > 1 while the analysis in [5] is restricted to d = 1). 
Compared to [5, Chapter 4], the main new result here is the observation that the 
contribution of these boundary conditions is given by a positive term (namely the 
first term in the right-hand side of (44)). 


Proof. The proof of this proposition relies on the divergence theorem applied to a 
well chosen vector field. Introduce the scalar function 

9 := X- Vx(t> + ydycf 


and the vector field 

= 9Vx,y(l). 

We are going to compute the integral of AiVx^y X by two different ways. The wanted 
identity (44) will be deduced by comparing the two results. 

Proposition 2.9 implies that Xx,y<i> £ (7^(0) and dnf = 0 on n 5H. This will 
allow us to justify all the computations done below. 


First computation. We want to exploit the fact that, since dn4> = 0 on R, one has 
X ■ n = 0 on R. To do so we begin by writing 



dnYx,y X dxdy 


X ■ ndS 


I an 


f X ■ n dS. 

Jan\R 


Since dQ \ R = {{x,y),x G Q,y = r/(x)}, by definition of G{r])’4}, the previous 
identity simplifies to 



diva;,^; X dxdy 


/ 9dn(pdS 

ldn\R 

f 9{x,r])^/TTWW dn4>\y=ydx 

>Q 



9{x, r])G{ri)'ip dx. 


Now, write 

Xxif = X x{(l>{x,r]{x))) = {X x(t)){x,ri{x)) + {dy(j)){x,ri{x))X xV- 

Since B = dy(j){x,r]), we get that {Xx4>)ix,r]{x)) = Xx'f’ — BXxXj- By definition of 
9, we deduce that 

9{x, 1]) = X ■ {Xxif - BXxV) + 'nB = X ■ Xxfj + {rj - X ■ Xx'n)B. 
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We thus end up with 


d\\x,yXdxdy= / {G{ri)'tf)){x ■ Vijj) dx + / {t] — x ■ Vx'n)BG{ri)'ip dx. (46) 

Jq Jq 


Second computation. Set 


W = IV 


x,y^ 


As can be verified by a direct computation, one has 

divx,y X = W + ^x- VxW + ^ydyW, 


and hence 


dwx,y X = div.T 


W 


W \ d-1 


x]+9y{^y] - 


W. 


Introduce the vector field 


Y = 


W 


2 \y^ 

Then the previous identity reads ddVx^yX = divx,yY — Consequently 

[[ diVa, y X dxdy = [ Y ■ ndS - 

Jan 2 

Now observe that, 


W dxdy. 


j^Y ■ ndS = ^ jjVx,y^f Q -ndS, 

[ Y-ndS = l [ W\y=. 

J d 


ldQ\R 
Therefore 


dx = ^ j {y-x- Vy)W\y=y dx. 


divx,y X dxdy = ^ j \'^x,y4'f -ndS 


\ f d \ 

+ 2 j^{'n-x-^d)y^\y=ri dx - — 


W dxdy. 


By combining this identity with (46) we obtain the wanted result (44). 


□ 


4 The main identity 


Theorem 4.1. Consider a solution {r],ip) G G^{[0,T]-,Hg{Q)xH^{Q)) of the Craig- 
Sulem-Zakharov system (20) with T > 0, s > d/2 + 2, d G {1)2} (see Defini¬ 
tion 2.W). Set 


0 = -ydtip - , 


and 


^{T) = L 1 L 2 J J Q{t,Li,X2)dx2 + ^ J Q{t,xi,L2) dx^ dt. 


1 
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Then the following identity holds 


^{T) = 'Ln + p + h + i2 + i3, 


(47) 


where Ti is the energy 


n = 


7 ; [ [V’G(r/)V^ + dx, 
^ JQ 


P is a positive integral (see Remark 3.2) given by 

rT 


P : = 


■=\ j [ (^\ndSdt 

^ Jo Jas^nan \yj 


and If are integrals denoting remainder terms: 
5 +2d 


h = 


h = - 


[ [ X,-h) dxdt, 

J JqxIO.T] 

{dy(j)){Vx'n ■ Vx<t>) dydxdt, 


8 JJqx[o,t] 
5 +2d r 


4 Jo JJn{t) 


d 1 


2 4 


I 3 = — [ -- / rj'ijjdx — r]{x ■ Vip) dx 


'Q 


t=T 


t=0 


'Q 


t=T 


t=0 


Before proving this result, let us deduce the following corollary. 

Corollary 4.2. Consider a solution (q^tp) E C^{^,T]]Hl{Q) x H^{Q)) of the 
Craig-Sulem-Zakharov system (20) for some T > 0 and s > d/2 + 2 with d E {1,2}. 
Assume that 

q> —/ r]dx = 0, (48) 

y Jq 

and that there exists two positive constants A, B such that 


B < 


T > 


5 +2d’ ^ - 2 - {5 + 2d)B 

(where max{Li,L 2 } = -hi */d = 1/ and 


^ ^ (2d + 3)max{Li,L2} ^ 

i H ;= A 

Vs 


Then 


sup ||V77(t)||^oo <-B, sup \\V'p{t)\\j ^2 < AVtH. 
te[o,r] iG[o,T] 


AS{T) > n. 


(49) 


(50) 


Proof. The proof is similar to the proof of Corollary 1.4. Firstly, notice that 

P> — If \Vx4>\'^{t,x,—h) dxdt. 

2 J Jqx[o,t] 

Therefore P + R >0 provided that /i/2 + ((5 + 2d)/8)q > 0, which holds true if 
rj > — 4 / 1/9 and d = 1, 2. 
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On the other hand 


,,, 5 +2d rT 

\h\ < —^—sup||Va:r/(t, Oli, 


/o J Jn(t) 


\dy4>\ \'^x4>\ dydxdt < 


5 +2d 


BTn. 


Using the Cauchy-Schwarz and Poincare inequalities, one has 

< max{Li,L2} \\r]\^2 ||Va:V'|lL2 , 


/ rj'ip dx 
Jq 


IQ 


TJX ■ Vip dx 


< max{Li,L2} ||7y|^2 ||Va,V’|lL2 , 


Since \\rj \\^2 < ^/2'H/g, we deduce from the assumption 50 that 

l-^sl < + 0 max{Li, L 2 }^^^. 

Therefore it follows from (47) that 

^{T) ^ + 0 max{Li, L 2 }-^A ) U, 

and hence SS{T) >T-Lhy assumption on T. 

Notation. We write simply 


J dx, J dy, j 


dt 


as shorthand notations for, respectively, 

rr]{t,x) 

IQ J-h 


r rl 

/ dx, / dy, 

Jq J-h Jo 


dt. 


(51) 

□ 


The proof of Theorem 4.1 will be made in three steps. First, we exploit the Pohozaev 
identity obtained in the previous section. 


Lemma 4.3. There holds 


{dtrj){x ■ V'ljj) dxdt = P — 


d- 1 


\^x,y(t>\^ dydxdt 
+ J J {x ■ Vrj — r]){dt'ilj + grj) dxdt. 

Proof. Since dtg = G{g)Tf, it follows from Proposition 3.1 that 

d- 1 


{dtr]){x ■ Vif) dxdt = P — 


1 


\'^x,y4>\ dydxdt 


(52) 


+ o {v-x-^v)iy +B- 2BG{g)Tf) dxdt. 
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To prove (52), it is thus sufficient to prove that 


+ 3"^ - 2BG{ri)il)) = -dti^ - grj. 


To obtain this identity, we use the fact that 


G{r])Tp + Vt? • VV' 

1 + |V7?|2 ’ 


V = Vijj — BVr]. 


Consequently, G{r])il! = B — V ■ Vrj and 


1 ,2 _ HVrr^^ + Girjy^ 

2^ 2 1 + |V7/|2 


\v‘^ + BV-Vg-]^B\ 


so (53) follows from the second equation of (20). 


( 53 ) 


(54) 

□ 


Next we integrate by parts and use the equations for (r/, ■0) to simplify the expres¬ 
sions. 

Notation 4.4. Let / : Q x [0,T] ^ M and d = 1,2. If d = 2 then we set 
r= ([0,Li] x{L 2 }) X ({Li}x [0,L2]), 

and use the notation 

I fLi I rL2 

— / f{t,xi,L2)dxi + — / f{t,Li,X2)dx2 
Jo X ‘2 Jo 

Similarly, when d = 1 and /: [0, Li] x [0,T] ^ 0, we write 


fdSdt = L 1 L 2 / 
rx[o,r] Jo 


dt. 


f dSdt = Li / f{t,Li)dt. 
rx[o,r] Jo 


With this notation, one has 


^(T) =11 0 dSdt. 

lrx[o,T] 


Lemma 4.5. There holds 


0 dSdt = P- 


d-1 


rx[o,T] 


\^x,y4>\ dydxdt 


— (1 -I- d) / / dxdt — 


2 + d 


— / r/(x • V'ljj) dx 


t=T 


t=0 


grf' dxdt (55) 


Proof. Notice that for any scalar function /, one has 


divxifx) dxdt = II f dSdt. 

rx[o,T] 
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Then, integrating by parts, we find that 


{x ■'Vrj — rj){dtip + gr]) dxdt = 


rx[o,T] 


grj^ dSdt — 


2 + d 


gg^dxdt 


rx[o,T] 


gdti) dSdt — (1 + d) / / gdti) dxdt 


+ {dtg){x ■'S/'ip) dxdt — / g{x dx 


t=T 


t=o 


We obtain the wanted result by combining the previous identity with (52). □ 

Remark 4.6. Consider the case when d = 1 and set m{t) = g{t,Li). Then, Since 
dxg{t,Li) = 0 and dxip{t,Li) = 0, it follows from (20) that 

dt'tp{t,Li) = -gg{t,Li) + ^{G{g)'ilj){t, Lif. 

Since G{g)ijj = dtg, we conclude that 

dtipit, Li) = -gm{t) + ]^m'{tf. (56) 

This shows that 0 = ^ \gm{t)‘^ — m(t)m'(t)^]. 


The next step consists in computing the right-hand side of (55). The wanted result 
(47) will be a direct consequence of the previous results and the following lemma. 

Lemma 4.7. Set 


A :— (1 + d) 


gdt'il) dxdt + 


2 + d 
2 


grj^ dxdt + 


d- 1 
2 


'^x,y4‘\ dydxdt. 


There holds 


A = - 


j J ['ipG{g)ij + gg"^] dx + ^ j g'^l^dx 

g\'^x4'\‘^it, X, —h) dxdt 


t=T 

t=o 


5 +2d 


+ III {dycl)){Vxg ■ Vj;(/)) dydxdt. 


(57) 


Proof. Let a be a parameter to be determined. We split A as A = Ai A 2 where 

d- 1 


Ai = a / / gdtif dxdt + 


^2 = (1 + d — a) / / gdtfj dxdt + 


dydxdt, 
2 + d f f ^^2 


Integrate by parts and use the equation dtg = G{g)'ip to obtain 


—a / / gdtf^ dxdt = —a / / 'il’G{g)if dxdt + a / gif dx 


J' 
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and 


so that 


d-l 


Ai = 


\'^x,y4>\‘^ dydxdt = 


'ijjG{r])il> dxdt, 


d-l 


— a 


dxdt + a / rjijjdx 


(58) 


On the other hand, directly from the eqnation for ip (see (20)) and the definition of 
B (see (28)), one has 


An = a - 


d 


1 


1 


grj'^ + {1 + d - a) I I rj ( + ^(1 + iVr/p)^^ ) . (59) 


In addition, it is easily verified (see (54) and (28)) that 


-^IVV’P + ^(1 + |Vr?|2)i?2 = 1 [(9,0)2 - |V0p - 2(a,0)(V0) • Vg 


y=n 


(60) 


The next calcnlations rely in a crucial way on the analysis done by Benjamin and 
Olver ([7]) of the conservation laws for water waves. This means that, to compute 
A 2 , we will use two elementary identities. Consider a scalar function u = u(x,y) 
and a vector field / = f{x,y) with values in ]R2 (the time is seen as a parameter and 
we skip it). Firstly, one has 



ri{x) 


dyu{x, y) dy + 


Q J—h 



—h) dx. 


(61) 


Similarly, 

/ f{x, 7]) ■ Vrjdx + 
Jo 


If 

JQ J-h 


divx / dydx = 


/, 


(9Qx [—/i,+oo))n90 


f -vdS. 


(62) 


Then 


= / / {dyU — diVa; /) dydx + / tt(x, —h) dx + 


f ■ V dS. (63) 


JJ J J (dQx[—h,+oc))r\dO. 

Recall that Ax^y4> = 0 and that dnp vanishes on n dH. (see Proposition 2.9). 
Therefore, by applying (63) with 

« = |(V)'-||V 0 P, f = -y{dyP)VP, 
we deduce from (60) that 


j V + i(l + |Vr/|2)R2^ = \ jj 

+ w / |V0|^(x,-/i)dx. 

J Q 
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Then, it follows from (59) and (58) that 




+ a J rjijjdx 

1 + d — a 


+ 


+ 


(1 + d — a)h 


'ipG{rf)il) dxdt + ( a — - 


[{dycp)"^ — |V0p] dydxdt 
|Vi?!)|^(t, X, —h) dxdt. 


grf dxdt 


We now chose 



1 

4’ 


Then the coefficients in front of ff iljG{r])ip dxdt and f f gr/^ dxdt coincide. On the 
other hand, notice that, since the energy % is conserved {d'H/dt = 0), one has 


\xl)G{g)'ip + gg^^^ dxdt 


['tljG{g)'tp + gg'^] dx. 


We thus find that 

A = — JJ ['ijjG{g)'ilj + gg'^] dxdt 



^ 5 +2d JJJ _ |y^|2j 

+ (5 +^2d)d JJ |V(/,|2(i^a:,-/i)dxdt. 


(64) 


Consequently, to obtain the wanted result (57), we need only to transform further 
the sum of last two terms in the right-hand side of (64). To do so, we use again (61) 
applied with (for some fixed t) 

u{x, y) = -{y - g{t, x)){dy(l)){t, x, yf. 

Then u{x,g{t, x)) = 0 and u{x, —h) = 0 and hence dyudy = 0. On the other 

hand 

dyU = -2{y - g)(j)y(j)yy - i(f)y)‘^, 

SO integrating on y E [—h^r]{x)] and then on x E [0,Li] we obtain, remembering 
that (pyy — 

^ = jj jj 

Since V(/> ■ u = 0 on [dQ x [—d, -|-oo)) O dO., by integrating by parts we infer that 


0 = - 


{y-v){dy'^= 


4>) • -|- 2 


xd ■ ^x4*t 
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so 


0 = - 


dy{iy-v)\'^x(p\‘^) + jj |Va;(^|^ + 2 jj (pyVa^r] ■ V^4> 

(|Va;0p -(pl) - j{h + r])\Vx4)\‘^ix,-h) dx + 2 jj (j)yVxr] ■ Vx4>- 


This proves that 


— \^ x4>\‘^) dydx + / h\Vx<p\‘^ix,—h) dx 


= - j r]\Vx4>\‘^{x,-h)dx+ 2 jj (pyVxr] ■ Vxfpdydx. (65) 
Plugging this identity in (64), we complete the proof of the lemma. □ 


By combining (55) and (57), we obtain the wanted result (47). This completes the 
proof of Theorem 4.1. 


5 Proof of Theorem 1.2 


In this section we prove Theorem 1.2. Let y G C'g"(T‘^) be even in xi and in X 2 - 
Consider an initial data {r]o,'ijjo) such that 


Vo ■= X{x) ^ 

|n| + |m|<A'' 


COS TT 


nxi \ 


cos TT 


mx2 

L 2 

mx2 


I ^ ^ 2 f nxi\ ( mx2\ 

V 0 -=X{x) 2^ O^^COS (^TT —j cos 


|n| + |m|<7V 


We furthermore assume that the mean value of r/o is 0, which in turn implies that 
the mean value of rj(t, ■) is 0 for all time. As explained below in Remark 5.3, one 
can consider more general initial data. 


Since these initial data are C°° and periodic functions defined for x G we can 
directly solve the Cauchy problem in by means of Theorem 2.3. Let (r/, ip) G 
C'°°([0,T] X T'^) be the unique solution of the system (20) with initial data {^ 0 , 1 ^ 0 ) 
(the fact that rj, ijj are C°° follows from the fact that on can propagate the regularity). 
As explained in Section 2, since the initial data are even in xi,X 2 , one then obtains 
solutions of the equations in the tank by restricting 


In view of Corollary 4.2, to prove Theorem 1.2 it remains only to prove the following 
result. 

Lemma 5.1. For any (3 > 1/2, there exist Kq,c, k such that, for any N G N*, if 

I -I- \o‘^ I ^ cN ^ 

I ^nm I ' I “nm | ’ j 

then the following properties hold with 

A = KqN^. 
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Firstly, the solution exists on [0,T(74)] with 


T{A) = 4 


^ ^ (2(i +3)max{Li,L2} ^ 
i H — A. 


and secondly one has r]{t,x) > —h/2 for all {t,x) G [0,T] x and 

sup ||Vr/(i)||^oc < —sup \\V'ip{t)\\j ^2 < AVm. 
tG[0,T{A)] b + M tG[0,T(yl)] 


To prove this lemma, we shall use an interpolation argument and the following result 
about the Cauchy problem on large time intervals. 

Theorem 5.2 (from Lannes [15]). Let si > 4. For all integer s > si there exist 
three positive constants depending only on s such that, if 

£0 := ll(%> < -5o, (66) 

then the solution (r/, V') given by Theorem 2.3 exists on a time interval [0,T] for 
some T > c*/eo and 


ll(r?(i),V^(t))||^ — ^1^0- 


Proof. Let us explain how to obtain this result from Theorem 4.16 in [15]. 


Set £ = Meq where M is some large enough real number to be determined later. 
Let (ry, if) be the solution given by Theorem 2.3 and introduce 


Then 


g‘^{t, x) = ir/(L x), if%t, x) = ^V'(L x). 


(67) 


dtp'' = G{£p'')'if^, 

+ gp" + - 2(1 + ^ 


( 68 ) 


The proof is based on two facts. Firstly, introduce the energy 


£'{t) = ||V# 


2 -L 

H2 + 


E 

|q|<s 


L2 


+ 




H' 


where 

G{£p')fj'' + eVp'^ ■ VV’^ 

1 + |eVry'^|2 

Then the hrst fact we need to know is a consequence of Theorem 4.16 in [15]. This 
result asserts that, if s > si and if f^(0) is small enough, say smaller than <5i, then 
there is a constant C* independent of £ and such that the solution exists on a time 
interval [0, T^] with > C'*/£ and satisfies 


sup £'{t) < 1. 
te[o,re] 
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To introduce the second fact we need to know, we begin by recalling the following 
technical ingredient: since 3 > (i/2 + 1, one has (see [15, Theorem 3.15]) 

11/^2 < C{\\£ri^\\Hi) \\^'4 ’"\\h^ ■ (69) 

Then, by combining this estimate with the usual nonlinear estimates in Sobolev 
spaces one obtains that, if 

< 1, (70) 

then ||ii^(t)||j |^2 is bounded uniformly in e € (0, Ij. Then, using the triangle inequal¬ 
ity and the product rule in Sobolev spaces, it is easily seen that 

Cl < £^{t) < C 2 ||(??^(t), VV’^(t))||^^+i , (71) 

for some absolute constants ci, C 2 . 


We are now in position to prove the wanted result. Fix M = C 2 /S 1 . Then 

We can always assume that 6 ijC 2 < 1, so that the condition (70) is satisfied at time 
t = 0. This allows us to use (71) which implies that 

£‘(0) < c, ||(tf(0),vv.‘(0))||^.,.< ftT = «i. 

As explained above, this implies that the solution exists on a time interval [0, T^] 
with > C^^je = (C*(ii/C' 2 )eo and satisfies sup^gjQ T^(t) < 1. Now observe that, 
by definition of if £^{t) < 1 then the condition (70) is satisfied. This allows 

us to use the bound (71) to deduce that 




1 1 C 2 

<—£ = —Meq = -^£q. 
Cl Cl Cidi 


(72) 


This gives the desired result where the parameters c*,(io,C'i in the statement of 
Theorem 5.2 are given by c* = C^ 5 i/C 2 , Sq = (^i/C '2 and Ci = C 2 /{ 6 iCi). □ 


We are now in position to prove Lemma 5.1. For some reason that will be clear 
below, we then fix s such that 


s G N, s > 4, 


s > 


/g + 1/2 

2;S - 1 ■ 


Once s is so fixed, we consider c*,(5o,C'i as given by the statement of Theorem 5.2. 
Recall that we consider initial data {'r]o,'tpo) such that 


Vo ■= Xix) 

|n|-|-|m|<W 


cos TT 


nxi\ { mx 2 \ 

r—— cos TT—— 


Li 7 


^0-=Xix) olmCOs(-K^'\ 

1^1 1 AT \ J- / 


COS TT 


|n| + |m|<A^ 


L2 J ’ 


mx 2 

L 2 


(73) 
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Directly from the Leibniz’ rule, we see that for any smooth functions ui, U 2 , there 
holds \\uiU 2 \\fja < ||tti||^ys.oo \\u 2 \\hs. Then, with 


Ho ■— ll(%)'^'*/’o)|lL2(Td)x_H'-i(T'*) > 


one has 


I I (% ) V V'O ) I I 1 ^ ^ ^ ^ 


< C 2 N^+^Ho, 


for some constant C 2 depending only on s and on ||x||,ys,cx). Set 


£0 := C2HoN^+^. 


Let a G (0,1] be a small positive number to be determined and assume that 


Ho < 


aSp 


(74) 


Then £0 < aSp/N^. Consequently, £9 < 5o and one can apply Theorem 5.2 to infer 
that the solution of system ( 20 ) with initial data (?? 0 )V’o) exists on a time 

interval [0,T*] with 

n> — > (75) 

£0 aoo 

and 

sup \\{'n,V'il)){t)\\^snd^^Hs-ind)<Cieo = CiC 2 HQN^^^. (76) 

te[o,r*] 


Now we use an interpolation inequality in Sobolev spaces: If u = Afii + (1 — A)(T 2 
with A G [0,1], then 

||'lt||/fo- < C*(cTi,(T 2) II'^II ll'*^ll_f/°'2 • 


This yields 

||VV'(t )||^2 < C(-l/2,s - 1) ||VV’(t)||^^!\ ||VV'(t)||j73i^^ with A(s) = 
Consequently, it follows from (76) that 

l|VVV)llL 2 <CMV'’Wify'<‘'||VVV)ll()!( with /!{>) = i£±tG, (77) 

and C{s) = C(-l/ 2 ,s - l)(CiC 2 )i-^(*). 

To conclude the proof if remains to estimate ||VV^||^_i and Hq in terms of the 
hamiltonian Ti := ^ J + ^G( 7 /)'i/’] dx. We claim that there exists an absolute 
constant K such that 

Ho<KVn. (78) 

Let us assume this claim and conclude the proof. Set 

Ko := KC{s)/V2 , A := KqN^. 
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It follows from (77) and (78) that 


l|V.^(i)||„ < KC(s)JV'’<‘)\/H with /3(s) = li±i^. 


By definition of s, one has /3(s) < /3. Therefore, by definition of A, 

||vv^(t)|li2 < KCis)N^Vn < AVm. 

Eventually, we chose a so small that the lifespan T* (see (75)) satisfies 



On the other hand, if the initial data are of the form (73) with < cN then 


Ho < c^lJ^{2N + 1)N-^, 


so that the condition (74) holds for k > s + |+/3 and c small enough. This completes 
the proof of Lemma 5.1 and hence the proof of Theorem 1.2. 


Therefore, to complete the proof of Theorem 1.2, it remains only to prove the claim 
(78). Notice that, by definition of Hq and H, one has 



so Hq < K\/TL will be a straightforward consequence of the claim || VV'H _ i < Ky/Ti 
at time t = 0 (with a different constant K). It is thus sufficient to prove the first 
inequality in (78). To do so, remember that (see Remark 1.5) y/2H is larger than 
the L^(ll(t))-norm of S/x,y(t>{t) (by definition (2)). Then, one infers a control of the 
Hx ' -norm of the trace (see [15, Prop. 3.12]) where the implicit constant 

depends on \\'r]{t)\\ff 3 - Now notice that one has a uniform control of ||T?(t)||// 3 , namely 
one has ||^(t)||j |^3 < 1 as a consequence of (76), the assumption s > 4 and the fact 
that one can assume without loss of generality that CiEq < 1. Therefore the Hx - 
norm of the trace V' is uniformly bounded by K^/H for some absolute constant K, 
which implies the desired result ||V'0|| _i < Ky/Ti. 

Remark 5.3. Let /3 > 1/2 and denote by s the smallest integer satisfying 



Given N gN and ^4,0 > 0, we denote by X]\f{A,c) the set of functions (r/o,V’o) £ 
R|+i(T‘^) X R|+i(T‘^) such that 



and satisfying 



Then the previous proof shows that, for any ^ > 0, there exists c > 0 such that, for 
any € N, if {rj, ipo) belongs to X]\f(A, c) then the conclusion of Theorem 1.2 holds. 
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A Hamiltonian proof 


The proof of Theorem 4.1, which was given in the previous sections, relies on several 
cancellations. To understand these cancellations, we will give in this appendix a 
(formal) proof of (16) which exploits directly the hamiltonian structure of the water- 
wave equations. ^ To simplify notations we consider the case d = 1 and assume that 
Li = 1 (then T = M/(2Z) and we denote by H^{T) the Sobolev space of 2-periodic 
functions). 

Consider a solution (ry,V’) S C^{[0,T]-,H^{T)xH^{T)) of the Craig-Sulem-Zakharov 
system (20) for some T > 0 and s > 5/2 (remember that the subscript e indicates 
that we consider even functions). Set 

A := [ [ [{dtr])x{dx'ilj) - (9fV’)a;(5x??)] dxdt. 

Jo Jo 

We are going to compute A in two different ways, and the wanted identity (16) will 
be deduced by comparing the two results. 


First computation. Here is the main new ingredient. In this step we compute A 
using the Hamiltonian structure of the water waves equation. This computation 
relies on the works by Zakharov [24], Craig-Sulem [12] and Lannes [14]. 


Following Zakharov [24] and Craig-Sulem [12], we begin by using the fact that the 
water-wave system (20) can be written as 


where^ 



dt'4^ 


m 

6 r] 


n 



+ dx. 


(79) 


The key point is the following: we have 

rT 


A = 


/ a(t} dt with a := 

Jo Jo 


'5% dn 

+ xr]^ — 
dtp dr} 


dx, 


which means that 


a = lim - ip + ex'ipx) — 'H{r}, tp)] + lim - [H{r] + exrix,^^) — 7^(f?, ip)] ■ 
£—^o e £—>-0 £ 


Since G{r]) is self-adjoint, one has immediately 


lim - [H{r],tp + exipx)-FL{'n,ip)] = [{G{r])'tp){xipx) dx. 

e^o e J 

^This proof could seem shorter than the one given in Sections 3-4 but this is not the case: 
indeed, it uses some computations done in Section 4, it is not self-contained (we use a shape 
derivative formula due to Lannes) and also we do not try to justify rigorously the computations. 

®By assumption, t] and pi are 2-periodic and even in x. This is why it is sufficient to consider 
integrals over [0,1]. In particular, the bilinear mapping {u,v) u{x)v{x) dx is a scalar product 

on Le(R/(2Z)). In (79), the derivatives of H are understood for this scalar product (instead of 
u{x)v{x) dx), and this is why (79) holds even if H is only one half of the energy. 
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On the other hand 


lirn — y \jj] + exrjxf'— dx = J rjxrjxdx 

so the only difficulty is to compute 


lim — 
£—>^0 2 s 


['4’G{r] + £xrix)'ip — dx. 


To do so, we use a formula due to Lannes (see [14, 15]) which allows us to compute 
the derivative of G{r])'il> with respect to t]. This formula gives 


G{r] + exr]^)^; = G{r])il) - eG{ri){Bxri^) - ed^iV xt^^) + 0{£^), 

where, as above, B = dy4>\y=rj and V = dx4>\y=rj. Using again the fact that G(r/) is 
self-adjoint, we obtain 


lim — [ ['il^G{ri + £xrix)i^ — V’G(r/)V'] dx 
£->■0 2e J 


[{G{ri)'il;){Bxr]x) + xr]x)] dx. 


By combining the previous results, we get 
“ = 5^ r]xr]xdx + j {G{r])'ilj){x'il>x) dx - ^ j [{G{r])'il>){Bxr]x) +'ipdxiVxrix)] dx. 
Since V{t, x) = 0 for x = 1 (by assumption) one has 

j i/jdxiVxtjx) dx = ^ J {xil)x)Vrixdx. 

Then, gathering the terms in a different way, we find that 

« = 9 y VX'qx dx + ^ j {G{r])'tp)x{'ipx - Brjx) dx + ^ J [{G{r])'ijj + Vr]x){xi;x)] dx. 


Remembering that 
we get 


'tpx - Brjx = V, G{t])iP+ Vr]x = B, 


a = g J gxgx 2 j dx + - J Bxipx dx. 

Recalling that m{t) = r]{t, 1) and computing the first term, we conclude that 


A = 


■ y m{ty 


dt — 


rf dxdt + 


1 


{G{r])^p)xV -|- Bxijjx dxdt. (80) 


Second computation. We now compute A using integration by parts and the equa¬ 
tions for rj and y. This second computation is not new. However, since it does not 
appear explicitly in the previous section, for the sake of readability, we redo the 
analysis. 
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Directly from the definition of A, by integration by parts, one has 

A = —J dt + j Tjxipxdx + JJ rjOti/j dxdt. 

It is convenient to split the last term as 

r]^t^p dxdt = - JJ rjdtip dxdt + - JJ rjdt'ilj dxdt. 

To handle the first term in the right-hand side, we integrate by parts in time and 
replace dtrj by G{r])ip. For the second term, we use the equation for ip written under 
the form (see (54)) 


dtip = -gr] - - EVd^r] + 


Then one gets 


Iff If t=T 

rjdtipdxdt = —- / / '4)G{rf)'ip dxdt + - rjipdx 


t=o 


35 


r]^ dxdt - r] -|- 2BVdxr] — dxdt. 


Recalling that dt'ip{t, 1) = —gm{t) -|- ^m'{tY (see (56)), we end up with 


A = y 'ygmftY — -rn{t)m!{tY^^ dt 

1 


^ YG{rj)ip dxdt —^ 11 dxdt 


+ 


1 / 


gif dx 


t=T 


t=o 


+ 


J xgd^ 


Ip dx 


t=T 


t=0 


g (r2 + 2BVdxg - B"^) dxdt. 


(81) 


Conclusion. By combining (80) and (81) we obtain that 

- f [gm{tY — m{t)m'[tY] ~ 7 f f Y’G{g)Y dxdt + ^ dxdt 




t=T f 
gY dx — xgdxY dx 
t=o ,/ 


t=T 


t=o 


+ R 


with 


(82) 


R := 


^g {V‘^ + 2BVdxg - B"^) + ]-{G{g)Y)xV + ]-B{xYx) 

o Z Z 


dxdt. 


Since 


B — dy(p\y=ri: — dx(p\y=rii G{jf)Y — gxdxY)\y='r)i 

and since Yx = {dxY + gxdyY)\y=y, one can write R as 

R= u{t,x,g{t,x)) dxdt + // f{t,x,g{t,x))gx{t,x) dxdt 
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where 


u{x,y) = + x(a^(/))(ay(/>), 

f{x,y) = ‘^y{d,;(l)){dy(i>) - + ]^x{dy(i>f. 

Then one has 

dyu - dxf = - {dy(t>f). 

Consequently, it follows from the identities (63) and (65) that 

^ = JJ 0 {dx(p)it,x,-Kf dxdt - ^ jjj {dx'n)idx(j)){dy(j)) dydxdt. 

The wanted identity (16) thus follows from (82). 
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